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ON THE FORMULATION OF THE AERODYNAMIC CHARACTERISTICS IN 


AIRCRAFT DYNAMICS* 

Murray Tobak and Lewis B. Schiff 
Ames Research Center 


1. INTRODUCTION 

As aircraft design has evolved toward more slender planforms in 
response to the demand for higher speeds, aerodynamic phenomena (e«g«, 
shock waves and concentrated vortex flows) have appeared that earlier 
were associated more with the flight of projectiles than with that of 
aircraft. The appearance of these flows has altered the aircraft’s 
aerodynamic behavior, and consequently its dynamical behavior, in ways 
that often are not predictable on the basis of UneoT'ized formulations 
of the aerodynamic force and moment system. It is generally agreed 
that consistently successful predictions of dynamical behavior under 
the new conditions hinge on incorporating an adequate nontinear depiction 
of the aerodynamic force and moment system within the equations of 
motion. 

A number of fundamental questions are raised by the necessity of 
replacing the linear aerodynamic formulation with a nonlinear one. 
Concepts such as stability derivatives, indicial functions, and super- 
position, which were the main working tools in the analysis and 
prediction of dynamical motions, are all based implicitly on a linearity 
assumption. Are these concepts at all viable when the linearity 

*Presented as a contribution to a course of lectures on "Aircraft 
Stability and Control" at the Von Karman Institute for Fluid Dynamics, 
Brussels, Belgium, May 12-16, 1975. 



assumption is removed, or must they be abandoned? How are time-history 
effects to be treated? For nonplanar motions, how are coupling effects, 
which could be neglected under the linearity assumption, to be treated? 
What about experiments? How are their requirements and their interpre- 
tation changed with the adoption of a nonlinear aerodynamic formulation? 

About ten years ago, the idea was introduced of using functional 
analysis as the vehicle for extending the linear formulation into the 
nonlinear domain [1,2]. The adoption of Volterra’s original conception 
of a functional [3] enabled the construction (without invoking a 
linearity assumption) of an integral form for the aerodynamic response 
that included the most general linear formulation as a special case. 

A framework was thereby established of sufficient breadth to enable 
answering some of the questions just posed. Initially, only planar 
motions were studied with a view toward answering the questions about 
the adaptability of the linear concepts and the treatment of time-history 
effects [1,2] . Later efforts focused on nonplanar motions, first of 
bodies of revolution [4-6], then of more general bodies [7,8], aimed at 
answering the question about coupling effects. Within the past few 
years, effort has been directed toward making the formulation more 
applicable to the treatment of dynamical motions of aircraft within the 
stall and post-stall regimes [9]. 

The occasion of this course provides the authors with a welcome 
opportunity to present a connected account of the theory’s main lines 
of development since its inception ten years ago. The paper begins with 
a brief review of the most general linear integral representation of the 
aerodynamic response for a planar motion involving two variables. An 
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examination of the principal failing of the linear integral forms leads 
naturally to the idea of replacing the indicial functions within the 
integrals by functionals, themselves dependent on the past motion. The 
formal (and rigorous) nature of this replacement is stressed to make 
clear that the treatment of more complicated motions involving any 
number of variables follows merely by an application of the same 
formalism. Then, simplifying approximations are introduced, aimed at 
reasonably specifying the degree to which the indicial responses can be 
cognizant of the past motion. The degree of cognizance allowed is shown 
largely to determine the generality of aerodynamic phenomena whose 
effects can be acknowledged within the scope of the resulting formulation. 
A simple specification leads to nonlinear forms of the aerodynamic 
response that can be recognized as practicable generalizations of the 
linear superposition and stability derivative formulations. Next, 
nonplanar motions are taken up and the question of coupling addressed. 
Simplification of the general formulation (allowing the indicial 
responses to have the same degree of cognizance of the past motion as 
before) yields a form for the total moment due to an arbitrary motion 
that can be compounded of the moment contributions from a limited 
number of well-defined characteristic motions. It is here that the 
particular importance of coning motion emerges as a motion characterizing 
the coupling problem. The characteristic motions in principle being 
reproducible in the wind tunnel, the result also provides the setting 
for a discussion of the kinds of wind-tunnel tests that, according to 
the formulation, may be required to enable consistently successful 
predictions of dynamical behavior. 
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COORDINATE SYSTEMS 


Three coordinate systems will be used. They have a common origin 
at the body's mass-center and a common axis x^ alined with a longi- 
tudinal axis of the body. 

Axes Xg, yg, Zg are body-fixed axes (Fig. 1). The flight 

velocity vector of magnitude V has components u_, v_, w resolved 

15 15 IS 

along x_, y_, z , respectively. Thus, 

B B ii 

'' - VV * "b^ 


y,q,M 



(a) CROSSFLOW PLANE (b) RESULTANT ANGLE -OF-ATTACK PLANE 

Fig. 1. Axes, angles, and velocity components In the crossflow and 
resultant angle-of-attack planes. 


Resultant angle of attack a Is defined by the flight velocity vector 

and the x„ axis. The plane formed by y„, z Is called the crossflow 
B B B 
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plane, illustrated in Fig, 1(a). The projection of a unit vector in the 
flight velocity direction onto the crossflow plane is a vector with 
magnitude 6; it will be called the (dimensionless) crossflow velocity 
vector. Reference to Fig. 1(b) gives 



The components of the body’s angular velocity relative to inertial space, 

resolved along x^, y^, z^, are p^, q^, r^, respectively. 

A second axis system x^, y, z is nonrolling with respect to 

Jd 

inertial space. Specifically, the x , y, z axes have an angular 

B 

velocity with respect to inertial space whose component resolved along 
the x^ axis is zero, while the components resolved along y, z are 
q, f , respectively. The angle ^ through which the body axes have 
rolled at any time t can be defined relative to the nonrolling axis 
system as 

* “ f 

•^0 

The angular Inclination X of the crossflow velocity vector 6 is 
measured relative to the nonrolling axis system, while ^ is the 
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angular inclination of the body axes from the crossflow velocity vector. 
With the aid of Fig. 1(a), the body roll rate is seen to be 

Pg = ^ ^ (4) 

The components p , q, r of the body's angular velocity vector resolved 
in the nonrolling axis system are related to those in the body axis 
system p^, q^, r^ through 

q + ir = e^^(qg + ir^) (5) 

Finally, an axis system x„, y, z will be called the aerodynamic 
axis system. Axis z lies in the crossflow plane and is alined with 
the direction of 6; axis y lies in the crossflow plane alined with 
a direction normal to the direction of 6. The components of the body's 
angular velocity resolved in the aerodynamic axis system p^, q, r are 
related to those resolved in the body axis system through 

q + ir = e^’^Cqg + ir^) (6) 


It will be noted (cf. Fig. 1(b)) that 6 and ip are no more than the 
polar coordinates of the dimensionless velocities w^/V and v^/V in 

the body axis system. Let w /V be called the angle-of-attack 

B 

parameter a and v^/V the angle-of-sideslip parameter 6; a and 6 
are related to the standard NASA definitions of angle of attack a and 
angle of sideslip 6 through 


tan a = w^/u^ = a/y 
sin S = Vg/V = e 
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and to 6 and ip through 


a + x3 = 6 


M 


( 8 ) 


The components of the aerodynamic force coefficient resolved along 

rs ^ 

the body axes x^, y^, Zg are Cy, C^, respectively, while the 

corresponding components of the aerodynamic moment coefficient (about 

^ A. A 

the mass-center) are C^, C^, C^, respectively. Analogously, the 

components of the force and moment coefficients resolved along the 

aerodynamic axes x^, y, z are C„, C^, and , C , C , 

o A 1 ^ i. m n 

respectively. Components of the aerodynamic moment coefficient resolved 
in the aerodynamic axis system are related to those resolved in the body 
axis system through 


h = 


(9) 


C + iC = e^’*'(C + iC ) 

m n ^ m n 


Corresponding relations between the components of the aerodynamic force 
coefficient are obtained by replacing m, n by X, Y, Z, respectively 
in Eq. (9) . 

To completely describe the state of a six-degree-of-freedom motion 

of a rigid body, it is necessary to specify the velocity and angular 

velocity vectors of the body. These may be expressed in terms of their 

scalar components resolved in the body-fixed axes u_, v_, w_, p_, q„, r_, 

D B 13 ij ij ij 

or equivalently by a, V, p , q , r„. Analogously, in the aerodynamic 

13 13 13 

axis system the motion is specified by the scalar variables 6, ij/, V, p_, 

B 

• • • 

q, r, or by V, X, q, r, since p is related to ip and X 

B 
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through Eq. (4). These are the principal variables on which the 
aerodynamic force and moment formulation must depend in general, 
although, as will be seen, certain simplifications are possible in special 
cases. 

3. DEVELOPMENT OF INTEGRAL FORMS 

For simplicity, the following (removable) conditions are imposed in 
all of the work to follow: (1) the aircraft is considered to be a rigid 

body which, prior to time zero, has been in steady flight for a long 
time; (2) at time zero, the aircraft begins a maneuver in which altitude 
changes are sufficiently small that atmospheric density and kinematic 
viscosity along the flight path remain essentially constant; (3) the 
aircraft's velocity along the flight path also is specified to remain 
constant. Hence, dynamic pressure, Mach number, and Reynolds number as 
measured along the flight path remain fixed throughout the motion. Under 
these restrictions, a general motion referred to the body axis system is 
described by the variables a, 3, p^,, r_; referred to the aerodynamic 

D D D 

% 

axis system, the motion is described by the variables 6, if;, X, q, r. 

The development of an integral form for the linear aerodynamic pitching- 
moment response to two of the variables will serve as a review of 
aerodynamic indicial functions and superposition. 

3.1 Linear Aerodynamic Pitching-Moment Response 

Let the aircraft begin a maneuver at time zero involving only two 
of the five flight variables on which the aerodynamic response depends 
(the others being fixed at zero). To fix ideas, the angle-of-attack 
parameter 6 and the pitching velocity q in the aerodynamic axis 
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system are chosen, so that the motion is planar. The extension to a 
more general motion involving all five motion variables will be obvious. 
For brevity, attention is focused solely on the aerodynamic pitching 
moment. All that is said, however, will hold as well for any other of 
the force or moment components that may exist merely on the substitution 
of that component for C^. 

Consider the aerodynamic pitching-moment response to the variations 
in 6 and q. It is usually admissible to break the variations into a 
large number of small step changes (cf. Fig. 2). In response to typical 


S(T) 


0 





A Cm 


AS 

1 

-(fiJ 




0 T t 0 T t 

Fig. 2. Summation of incremental responses. 


step changes A6 and A(qZ-/V) at time t, there is an incremental 

change in pitching moment AC ; it is measured at a fixed time t 

m 

subsequent to x. The assumption of linearity is now invoked, having, 
in the context of this report, the following meaning: AC^ is said to 
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be independent of 6 (t), q(x), and the past values of these variables. 

This enables writing AC in the form 

m 


M ~A(qi/V) ' 

The form implies that AC /A6 and AC /A(qZ/V) are derivable from 

m m 

linear differential equations. That they depend on elapsed time t t, 
rather than on t and t separately implies that the coefficients of 
the differential equations are independent of time. It is important to 
note the following distinction: the significance of the linearity 

assumption does not rest on the assertion that AC^ is linearly 
dependent on A6 and A(qZ/V), or that the contributions to AC^ from 
the two increments are linearly additive. Both of these assertions are 
always justified whenever A6 and A(q^/V) are so small that terms of 
0((A6)^, (AqZ./V)^, (A6AqZ/V)) can be neglected. The significance of the 
assumption lies in saying that the two ratios are independent of 6 
and q. Thus, no matter how large the values of S and q at the 
origins of the steps, the response functions AC^/A6 and AC^/A(qZ/V) 
are said to be the same functions, dependent only on t - t. The limits 
of these functions as A6 and A(q^/V) approach zero 


lim 
A 6’^0 


AC (t - t) 
m 

A6 




(t - t) 


AC (t - t) 

lim s - r* - r'i 

Hqim-H) A(qZ/V) 



are called the indlcial pitching-moment responses per unit step changes 
in 6 and q^/V, respectively. Every pair of steps in 6 and qt/V 
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beginning at a value of t less than t has a corresponding incremental 

pitching-moment response of the form Eq. (10) that contributes to the 

pitching moment at time t. The summation of these incremental responses 

to the steps that occur over the range of r from zero to t then gives 

the aerodynamic pitching-moment response at time t to the variations 
in 6 and q* As the indicial responses depend only on the time 
difference t - x, in the limit the summations take the form of the 
familiar convolution integral. The sum of the two contributions and the 
initial value of then give the total pitching moment at time t 


C (t) = C (0) 
m m 


/' 




t) 6(T)dT 


- r 

V Jn 


m 


(t - 


t) q(T)dx 


( 12 ) 


For constant V and within the assumption of linearity, Eq. (12) is 
exact. Since all values of 6 and q figure within the limits of the 
integrals, at time t depends on the whole past of 6 and q. 

3.2 Definition of a Functional 

The description of C (t) as a function that depends on all of the 

m 

past values of the argument functions 6 and q corresponds mathemati- 
cally to Volterra*s description of a funationat [3] , If Volterra’s 
square-bracket notation is adopted, the assertion that ^ 

functional is indicated thus, 

C^(t) = G[6(5),q(0] (13) 

where it is understood that 5 is a running variable in time, ranging 
over the interval zero to t. 
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In brief, just as an ordinary function f(x) assigns a number to 

each X for which it is defined, a functional F[y(5)] assigns a 

number to each function y(?) of the set of functions (all of which 

are defined in some interval a < 5 < b) for which the functional is 

defined. Thus, Eq. (13) may be interpreted as follows: given any pair 

of functions 6(g),q(C) out of the collection of all such pairs defined 

in the interval 0 1 5 < the functional G assigns a number to 

C (t). It is the idea of a functional, with its formal recognition of 
m 

the influence of past events on present behavior, that provides the key 
to generalizing the linear integral form. 

3.3 Nonlinear Aerodynamic Pitching-Moment Response 

An heuristic argument will show how the most general linear form 

for the aerodynamic response ^jj^(t) to two of the flight variables 

(Eq. 12) can be freed of the linearity assumption. A more rigorous 

mathematical development yielding the same result is available in Ref. 1. 

Attention is directed to the incremental form for (10) > 

where the principal limitation imposed by the linearity assumption first 

appears. As already noted, this limitation does not hinge on the 

assertion that AC (t) is linearly dependent on A6 and A(q^/V), or 
m 

that the two incremental contributions are linearly additive. Both of 
these assertions remain in force in the more general development. The 
limitation is simply that the indicial responses are said to be independ- 
ent of the past values of 6 and q. Now, as described more precisely 
below, an indicial response to a step change in a motion is formed by 
taking a difference between the responses to two motions whose histories 
differ only by the step imposed on the second motion. Hence, the motions 
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prior to the step are identical. When the difference is taken in the 
linear case, the influence of this identical past cancels identically. 

In the absence of linearity, exact cancellation of the past influence is 
not to be expected, so that some remnant of that influence must be 
evidenced in the behavior of the indicial response. That is to say, the 
indicial response must be a functional. The way to eliminate the 
limitation of the linearity assumption in Eq. (10) thus becomes clear: 
replace the indicial functions by functionals, themselves dependent on 
values of 6 and q prior to the origin of the steps. That this 
replacement can be rigorously justified is, in fact, the principal result 
of the mathematical development in Ref. 1. 

If the notation already introduced for a functional is followed, the 
designation of the indicial responses as functionals is indicated thus, 


AC (t,x) 

lim — r6(0.q(5);t,T] 

A6->0 ° 


AC (t,T) 

■ S'''"-'’'”'''"' 



where, as before, t is the time at which the increment AC is 

m 

measured, x is the time at which the steps originate, and 5 is a 
running variable in time over the interval zero to x; that is, over the 
past for steps beginning at ^ - x. 

Although the replacement of the indicial functions by functionals 
appears to be largely a matter of notation, the change is far-reaching 
and, for example, requires a more precise description of the formation 
of indicial responses than was necessary in the linear case. Two motions 
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AS 




Fig. 3. Formation of indicial responses. 
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have to be considered (cf. Fig. 3): first, beginning at 5 = 0, the 

aircraft is made to execute the motion under study 5(5), ^ . 

At a certain time t, the motion is constrained such that the values 
of the flight variables existent at time x, that is, 6 (x) , y — 
remain constant thereafter. The pitching moment corresponding to this 
maneuver is measured at a time t, subsequent to x. Second, the aircraft 
is made to execute precisely the same motion, beginning at 5=0 and 
constrained in the same way at 5 = x, except that at the latter time, 
one of the variables 6 or qZ/V is given an incremental step A6 or 
A(qZ-/V) over its value at 5 = t. Hence, if it is 6 that is given an 
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increment A6, the values of the flight variables for all times 
subsequent to x are 6(x) + A5, pitching moment 

corresponding to this maneuver is again measured at time t. The 
difference between the two measurements is divided by the incremental 
step A6 or A(ql/V); the limit of this ratio as the magnitude of 
the step approaches zero is called the indicial pitching-moment response 
at time t per unit step at time x of one of the two flight variables 
6 or qZ^/V, As indicated in Fig. 3, since the two motions prior to 
Z - X are identical, the ratio must be identically zero for all 5 < x. 
At 5 = T, a discontinuity in the ratio is permissible, reflecting the 
discontinuous change in one of the motion variables. For all values 
of 5 > X, the ratio must be continuous. 

Formed as described above and as defined in Eq. (14), the indicial 
responses are suitably generalized to be free of dependence on a linear- 
ity assumption. With Eq. (14) replacing Eq. (11) in Eq. (10), the 

summation of incremental responses to yield an integral form for C (t) 

m 

follows as before. The resulting generalization of Eq. (12) is 

^ Cm^t6(5),q(5);t,T] -^6(x)dT 

*'o 

+ C^^[6(C),q(C);t,x] q(x)dx (15) 

Equation (15) , applicable to a planar motion involving arbitrary varia- 
tions in 6 and q, is believed to be of sufficient generality to 
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provide a framework for the study of a wide variety of nonlinear 
aerodynamic problems. Derivation of the integral forms corresponding 
to motions involving more than two variables follows simply by an 
application of the same formalism. In the next section, simplifying 
approximations are introduced into Eq. (15) which make the form more 
practicable at the price of narrowing the classes of aerodynamic behavior 
that can be acknowledged within its scope. 

4. LIMITATIONS AND SIMPLIFICATION OF THE NONLINEAR INTEGRAL , FORM 

Conditions that flight-path properties remain constant during the 
motion and that the aircraft be a rigid body have been imposed at the 
outset. Imposing constant flight-path properties excludes from 
consideration the influence on motions of very large accelerations or 
variations in atmospheric density, such as might occur, for example, 
during atmospheric reentry. Omitting structural variables rules out the 
possibility of treating the buffeting problem, which involves interactions 
between the elastic airframe and random aerodynamic fluctuations [10]. 
However, the presence of fluctuations themselves can be acknowledged 
within the framework already established by adopting the notion of 
ensemble averaging (cf. Ref. 9). The remaining assumptions are of two 
main classes: first, fundamental assumptions associated with the use of 

functional analysis to develop the general integral form for the 
aerodynamic response; second, simplifying assumptions associated with 
the reduction of the general integral form to more practicable forms. 
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4 . 1 Fundamental Assumptions 

Despite their generality, the nonlinear indicial responses defined 
by Eq, (14), and hence the corresponding integral form for the aero- 
dynamic pitching moment, Eq, (15), already contain implicit assumptions 
that limit their applicability. These are principally that the indicial 
responses must exist for all values of their arguments (for 5 > ^i^d 
must be unique. Evidence for the existence of indicial responses is 
necessarily inferential, since it is doubtful whether they will ever be 
measured directly. However, the nonexistence of indicial responses 
follows from the nonexistence of their steady-state values, which are, 
of course, more amenable to measurement. In this respect, assuming the 
existence of the indicial responses for all values of their arguments 
(with 5 > t) strictly speaking implies the exclusion of cases where the 
variation of steady-state pitching moment with 6 (e.g., at q 0) 

becomes discontinuous either in its magnitude or slope at certain 

ACm(«) 

isolated values of 6, since lim ^ — will not exist at these 


A6-K) 


points. Such cases are known to characterize certain types of stall 
behavior (cf., for example, Ref. 11). Although these cases can be 
treated by an appropriate addition of jump conditions at the isolated 
points, for simplicity they will be excluded from further consideration 
here. Assuming uniqueness implies the exclusion of cases where more 
than one steady-state response to the same maneuver is possible; more 
precisely, cases must be excluded where the steady-state responses 
to repetitions of the same maneuver can be multivalued and probabilistic. 
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An example of aerodynamic behavior that must be excluded on this basis 
is the steady asymmetric vortex flow above a body of revolution at 
large 6 when it is a matter of chance whether the vortices form a 
left-hand or a right-hand pattern. Given these fundamental limitations 
on the applicability of the general integral form for Eq. (15), 

simplifying assumptions still need to be attached in order to reduce it 
to more usable forms. 

4.2 Simplifying Assumptions 

In the form Eq. (15) the indicial responses within the integrals 
are themselves functionals, depending in general on the whole past 
history of the motion 6(5), q(C). This makes the further use of the 
form exceedingly difficult, since the history of the motion normally is 
not known in advance but rather is desired as the solution of the 
equations of motion. Thus, when the past history is unspecified, the 
functionals also are unknown beforehand. Simplification of Eq. (15) 
hinges on replacing the functionals by appropriate functions whose 
dependence on the past is denoted by a limited number of parameters 
rather than by continuous functions. If 6(5),q(5) can be considered 
analytic functions in a neighborhood of 5 = t (corresponding to the 
most recent past for an indicial response with origin at 5 = t) , in 
principle their histories can be reconstructed from a knowledge of all 
of the coefficients of their Taylor series expansions about 5 = 

Thus, since 6(5),q(5) are equally represented by the coefficients of 
their expansions, the functional, with its dependence on 6(^),q(5), 
can be replaced without approximation by a function with a dependence 
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on all of the coefficients of the expansions of 6(5), q(?) about 5 = t. 
The indicial response example, can be expressed as 

= Cnig(t,T;6(x),(S(T) q(T),q(t). . (16) 

(The reversal in the order of the dependencies anticipates the diminished 
role played by 6 and q in the succeeding analysis.) Now physical 
reasoning suggests that the indicial response should have "forgotten" 
long-past events and thus should depend mainly on events in the most 
recent past. If this is assumed to be true, then so far as the effect 
of the past on the indicial response is concerned, the form of the past 
motion just prior to the origin of the step might just as well have 
existed for all earlier times. Hence, at most only the first few 
coefficients of the expansions of 6(5),q(5) need be retained to charac- 
terize correctly the most recent past, which is all the indicial response 
remembers. Retaining the first two coefficients of 6(5), for example, 
implies matching the true past history of 6 in magnitude and slope at 
the origin of the step, thereby approximating 6(5) by a linear function 
of time 6(5) 6(t) - 6(t)(t - 5)* With an approximation of this order 

in force for both 6(5) and q(5) in the indicial responses, the integral 
form replacing Eq. (15) becomes 

^ q(x)dT (17) 

This form, while considerably more tractable than Eq. (15), is still 
sufficiently general to allow the treatment of motions involving 
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hysteresis effects. Retaining a dependence on 6(x), for example, allows 
assigning different indicial responses to a step at a single pair of 
values 6(T),q(x), depending on the magnitude and sign of 6(t). It is 
permissible, for example, to distinguish between indicial responses 
where S was increasing or decreasing prior to the step. This will be 
valid when physical reasoning suggests that the particular response 
attached to a step is a deterministic result of the past history rather 
than the probabilistic result of an interaction with a random fluctuation. 
Although implications of the use of Eq. (17) will not be explored further 
in this study, the equation is believed to be both tractable enough and 
of sufficient scope to provide a framework for the study of rapidly 
varying maneuvers (for example, the rapid pull-up), where hysteresis 
effects governed by rate-dependent flow phenomena are known to be 
present in the aerodynamic response. 

When deterministic hysteresis effects are absent, and if, addition- 
ally, the assumption of a slowly varying motion is introduced, then the 
dependence of the indicial responses on 6(T),q(T) will not be signifi- 
cant (cf. Fig. 4). So far as the indicial responses are concerned, 
omitting these dependencies in Eq. (17) implies that the motion prior to 
the origin of the steps is being approximated by the time- invariant 




Fig. 4. Slowly varying motions. 


20 



motion 6(?)«5 5(x),q(5) « q(T). The indicial responses at any value 
of elapsed time, now dependent only on the magnitudes of 6 and q 
just prior to the steps, must not only be continuous functions of 

but henceforward also singte-valued functions of 5(T),q(x). 
Further, with a given time- invariant past motion and with the already 
assumed constant flight-path properties, clearly an indicial response 
must have the same value after a given time has elapsed subsequent to 
the origin of a step no matter when the step occurs. That is, just as 
in the linear case, the indicial response must be a function of elapsed 
time t - X rather than of t and x separately. Finally, then, for 
slowly varying motions and within the additional restriction on the 
indicial responses of single-valuedness with respect to 6(x),q(x), and 
t - X > 0, a much more specific form of Eq. (15) may be written, still 
capable of embracing a fairly broad range of nonlinear aerodynamic 
problems. It is 

+ Y A q(T)dT (18) 

Although the form of Eq. (18) represents a great simplification over 
that of Eq. (15), the equation still includes the full linear form 
(Eq. (12)) as a special case. Equation (18) is the basic integral form 
underlying the subsequent simplified formulations. 
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5. APPLICATION OF THE SIMPLIFIED NONLINEAR INTEGRAL FORM TO PLANAR 

DYNAMIC STABILITY STUDIES 

Equation (18) is now applied to the study of aircraft dynamic 
stability. The rigid-body motions of aircraft are normally oscillatory, 
and moreover, the oscillations are generally of very low frequency. 

Several analytical benefits accrue from the latter fact. First, since 
the motions are slowly varying, the assumptions underlying Eq. (18) are 
particularly well grounded in this application. Second, Eq. (18) can 
be further simplified. The simplification, which in effect reduces 
Eq. (18) to an equation correct to the first order in frequency, parallels 
that realized in the linear case in the application of Eq. (12) to 
stability studies [12]. 

5.1 Stability Coefficient Formulation 

The reduction of Eq. (18) to first order in frequency is carried 
out in detail so that, later, the analogous reduction of the more lengthy 
equations for nonplanar motions will need only to be Indicated in passing. 
Equation (18) is first rearranged to give a more convenient form. From 
physical considerations, the indlcial responses must approach steady- 
state values with Increasing values of the argument t - x. To indicate 
this, the following notation is introduced (the notation parallels that 
of Ref. 12) : 


Cm5(<^ " T;6(x),q(T)) 

Cmgl 

^«;6(x),q(x)^ 

- ?Jt - x;6(x),q(x)) ] 





Cmq(t - T5'S(x),q(x)) 

^qi 

^“;6(x) ,q(x)^ 

- Fi^^t - x;6(x),q(x)j j 
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where 


Cnig(«;6(T),q(T)) 


Cm^^;6(x),q(T^ 


rate of change with 6 of the pitching-moment 
coefficient that would be measured in a steady 
flow, evaluated at the instantaneous value of 
6 (t) with q fixed at the instantaneous value 
q(T) , 

rate of change with q of the pitching-moment 
coefficient that would be measured in a steady 
flow, evaluated at the instantaneous value q(x) 
with 6 fixed at the instantaneous value 6 (t) . 


The functions and are termed deficiency functions; they tend 

to vanish with increasing values of the argument t - t. When Eqs. ( 19 ) 
are inserted in Eq. ( 18 ), the terms involving the steady-state parameters 
form a perfect differential which can be immediately integrated - 
Equat ion (18 ) b ecomes 

C^(t) = C^^~;6(t),q(t)) - J' Fg^t - x;6(T),q(T)j ^ 6(x)dx 

- x;6(x),q(x)) q(x)dx (20) 

total pitching-moment coefficient that would be 
measured in a steady flow with 6 fixed at the 
instantaneous value 6(t) and q fixed at the 
instantaneous value q(t). 
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where 

C^(~;6(t),q(t)j 


Equation (20) is a form of Eq. (18) particularly amenable to 
approximation. Let it be assumed for illustration that the angle-of- 
attack parameter 6 is essentially a harmonic function about a constant 
mean while the angle of pitch 0 (0 E q) is a purely harmonic function. 
That is 


6 


6 + 
m ^ 


io)t 

e 


Q ^ Q 


0 


io)t 

e 



. 

Then clearly, since q = 6, q itself will be of first order in 
frequency to. Hence q will be small for all values of time, and powers 
of q higher than the first will be of second and higher orders in fre- 
quency. Therefore, for any given values of t or x, it is permissible 
to expand the terms in Eq. (20) in a Taylor series about q = 0 and to 

discard terms containing powers of q higher than the first. Terms in 

. • 

q and 6q likewise may be discarded as they will be of second order in 
frequency. The result of the expansion is 


- Fg^t - t;6(t),o)^ 6(x)dT 

•'ri 


( 22 ) 


Definitions of C^^“;6(t),oj and (t) ,oj follow from those given 

earlier with the substitution of q(t) = 0. The first two terms are 
clearly the nonlinear counterparts of the terms SCmg (") and (qZ-/V)Cjjj^(«>) 
that appear in linear analyses based on the stability derivative concept. 
Therefore, the integral, when also reduced to the first order in frequency. 
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is anticipated to be the nonlinear counterpart of the term • 

This reduction is taken up next. 

With the change in variable t - t = the integral becomes 

rt / \ iw(t-Tj) 

1=1 F 3 ^T^; 6 (t - Tj),Ojio) 6 Q e dx^ (23) 

•'o 

which may be rewritten 

. -icox , 

F3(xi;6(t - x^),oje dx^ (24) 

Practically speaking, the deficiency function essentially vanishes after 
a relatively short period of time has elapsed. Let the value of x^ at 
which F essentially vanishes be t , and consider events at a time t 

o Si 

sufficiently removed from the start of the motion that t > t^. Then 
the upper limit in Eq. (24) may be replaced by t , whereupon, with x 
bounded and (a) small, the harmonic function may be expanded in powers 
of w. Since 6 is itself of first order in oj, however, only the first 
term in the expansion, unity, contributes within the order of the 
approximation. Moreover, with respect to the parameter 6 (t - a 

further simplification can be realized when the condition t > t^ is 
invoked, for then 6 (t - « 5 (t) . The integral reduces to 


where, as anticipated, the integral term within parentheses may be 
identified with “Cmg • 
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Hence, just as in the linear case [12], to the first order in 
frequency 0^^ is proportional to the area of the deficiency function, 
now, however, evaluated at and dependent on the particular value of the 
angle-of-attack parameter 6(t) under consideration. The nonlinear 
counterpart of the linear stability derivative formulation thus becomes 

^m^^^ “ Cmq(«>;6(t),o) + 6 (t) (fi (t)) (26) 

where 

Cm^(^(t)) =-f / F3(x^;6(t),0)dT^ (27) 

•'o 

5,2 Interpretation of C 2 ^ 

In Eqs. (26) and (27), just as in the linear theory, , the term 
that accounts for the past, appears as an integral of the deficiency 
function, A physical argument will show why the term, or more directly, 
the analogous normal-force coefficient continues to appear in this 

way. In Fig. 5, the aircraft on the left has been sinking, without 
pitching, for a long time at a constant rate. A force P must be 



Fig. 5. Interpretation of • 
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applied to maintain the constant rate. The work done by the applied 

force over an arbitrarily large time interval zero to t is 

a 


Wk^ 



[W + Z(oo;6Q>]dT 


( 28 ) 


where W is the weight of the aircraft and -Z(«>;6 q) is the steady-state 
normal force due to the constant angle of attack. Now, as shown on the 
right side of Fig. 5, let the same aircraft experience a step change in 
6 at time zero and then undergo the same motion as in the first case. 

The work done over the same time interval by the force applied to 
maintain a constant rate is 


»/ 


Wk 2 = V6o I [W + Z(x;6o)]dT 

*0 


(29) 


The difference in work done is 


r 


Wk^ - Wk^ = V6 q I [Z(oo;6q) - Z(T;6Q)]dx 


(30) 


After identifying t^ with the time required for Z(t;6q) to reach 
steady state, one sees that the integral is the area enclosed by the 
indicial normal-force response curve and its steady-state value. That 
is to say, it is the area of the normal-force deficiency function, and 
it is therefore proportional to • The energy of the aircraft is 

the same in both cases, since it undergoes the same motion. The energy 
expended by the applied force is different in the two cases. The 
balance of energy, which is Eq. (30), therefore must have been given to 
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or taken from the fluid. The term is a measure of the energy given 

to or taken from the fluid whenever the angle of attack changes from one 
level to another. Because this assertion holds regardless of the 
magnitude of the angle of attack, it is reasonable that (and 

should continue to appear as an integral of a deficiency function, even 
in the nonlinear analysis. 


5.3 Theoretical Evaluation of Stability Coefficients 

Of the three coefficients in Eq. (26), the first is the familiar 
static pitching-moment coefficient due to angle of attack. Since it is 
of fundamental importance in aircraft design, it has been studied 
extensively; a large body of both theoretical and experimental results 
may be called on to define it in particular cases. No such body of 
information exists for the remaining two coefficients. The specific 
problems which must be solved for these two coefficients are character- 
ized by the motions illustrated in Fig. 6. Just as for the first term, 
solutions for may be derived from a steady-state equation. 



S =CONSTANT 
^ = CONSTANT— 0 



(a) MOTION FOR Cm-(S) (b) MOTION FOR Cmg(®:S,0) 

8 ^ 


Fig. 6. Definitive motions for calculation of stability coefficients. 
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In fact, the theoretical problems for and for the static pitching- 

moment coefficient due to angle of attack differ little in principle; 
the computational methods that have been devised for the numerical 
solution of the latter problem should be largely transferable to the 
solution of the former. Solutions for other hand, must be 

derived from a time-dependent equation. Except for those cases where 
one of the three space coordinates can be omitted or where a similarity 
condition holds (e.g., two-dimensional flow, conical flow), the necessity 
of treating an essentially four-dimensional problem apparently places it 
beyond the capacities of even present-day automatic computing machines. 

At any rate, the authors are unaware of any attempts to treat the problem 
for with an allowance for the nonlinear dependence on <S . In view 

of this, some of the physical features of the indicial response and of 
the integral form for will be brought out that might assist in 

making at least order-of -magnitude estimates of • 

The two-dimensional wing is chosen for demonstration purposes since 
the important features of sound-wave propagation are most easily illus- 
trated in this case. A subsonic flight speed is specified because here 
many of the results from the linear theory do not require modification 
in principle and thus can be adapted to guide the drawing of a physical 
picture of the indicial response. By means of these adaptations, the 
results may have more general bearing. Consider the boundary conditions 
corresponding to the motion illustrated in Fig. 6(a), where the flight 
speed V is subsonic. Let the wing move away from a coordinate system 
that is fixed in space at the position of the mass-center at time - 0 
where the step change in the boundary conditions occurs (cf. Fig. 7). 
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Fig. 7. Boundary conditions for indicial loading. 

The distance traveled by the mass-center along the flight path is 
measured by a coordinate s^. Since the mass-center moves at constant 
speed, the trace of its path plotted against time is a straight 

line. This is shown in Fig. 7, where the maximum projections of the 
leading and trailing edges are also shown, parallel to the trace of the 
mass-center. For x^ < angle-of-attack parameter 6 is constant 

at 6(t). A step change in 6 occurs at x^^ = 0, so that for x^ > 0, 
the angle-of-attack parameter is 6(t) + A6. Due to the impulsive 
change in 6, the loading on the wing at x^ = 0 also undergoes a 
sudden change. The physical situation at this instant corresponds 
essentially to that described by piston theory (cf., for example, 

Ref. 13), which should give a reasonably accurate estimate of the initial 
change in loading. Reflecting the local character of the instantaneous 
response, the initial change in loading is essentially constant across 
the chord, so that the corresponding center of pressure is located very 
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near the mid-chord. Also at = 0, the sudden change in flow conditions 
causes disturbances to be propagated from every point on the chord line. 

Of particular Importance are those from the leading and trailing edges. 
Each disturbance is propagated at the local speed of sound, so that, on 
a plot such as shown in Fig. 7, the zone of its influence is bounded by 
projections of an approximately conic surface whose origin is the point 
of the disturbance. As shown on Fig. 7, traces of the waves from the 
edges divide the wing into a number of distinct regions. Points in 
region (1) have not yet been made aware of the changed conditions at the 
leading and trailing edges by the arrival of the sound waves, and hence 
the loading in this region remains essentially unchanged from that 
existent at = 0 (cf. Fig, 8). This loading gradually disappears as 

the propagation of the two sound waves announces the new conditions to 
increasing portions of the wing, and it disappears completely at 
x^ = t^. Points in region (2) have been made aware of new conditions 



Fig. 8. Indicial loading on two-dimensional wing at subsonic speed. 
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at the leading edge but are as yet unaware of changed conditions at the 
trailing edge. The reverse is true for points in region (3). Points 
in region (4) have become aware of changes at both the leading and 
trailing edges, whereas points in region (5) in addition have become 
aware of the reflection from the trailing edge of the initial wave from 
the leading edge. Not long after the first reflection, the loading on 
the wing has essentially assumed the form of its ultimate steady-state 
loading (Fig. 8) and only its amplitude increases thereafter with 
increasing tj, asymptotically approaching the steady-state amplitude. 
Therefore, as shown in Fig. 8, the center of pressure has essentially 
reached its steady-state position while the loading is still increasing 
in amplitude. The above behavior suggests breaking the indie ial 
pitching-moment function into two separate contributions as shown in 
Fig. 9. The first variation represents the pitching-moment contribution 
of the integrated loading in region (1). Accordingly, it begins with 
the value ^0; 6 (t) , oj and vanishes at = t^. The second variation 

reflects the lumped contributions of the integrated loadings in the 
remaining regions; in conformity with the loading, its initial value 


Crr, 


8 


Crr^ 


8 


(t|: S(t), o) 



Fig. 9. Breakdown of indlcial response. 
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is zero while its end value, essentially attained at x, = t , is 
(t) ,o). The sum of the two contributions is the indicial 
pitching-moment function ^x ;6 (t) ,0 The end values of the two 
contributions are evidenced explicitly by introducing the normalized 
functions f^ and f^ as shown in the figure. This enables writing 
the indicial response in the form 

Cmg(Ti;6(t),o) = Ci„^(0;6(t),0)f^(T^;6(t)) 

+ (31) 

where f^ and f^ vary within the limits zero and unity. Then the 
deficiency function is 

F3(T^;6(t),o) = Cni5(~;6(t),o)[l - f2(T^;6(t))] 

- CV^(0;6(t),0)fi(Ti;6(t)) (32) 

SO that the integral for takes the form 

I f^(x^;6(t))dT^ 

0 

- Cm5(“;'5(t),o) I Fl - f2(x^;6(t))JdT^ 

Jo 

(33) 

Equation (33) clearly shows the strong dependence of on the 

steady-state parameter Cm5(”;5 (t),0^. Further, the sign relation 
evidenced in Eq. (33) should be noted. Since the integral involving f^ 
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normally will be positive, the equation indicates that as (t) ,0j 

becomes more negative (i.e., statically more stable), becomes more 

positive (i.e., dynamically more unstable). 

Since and are normalized functions, cases are anticipated 

where their dependence on 6 will not be significant. For example, the 
process of arriving at the form Eq. (33) suggests that it might be 
reasonably accurate to use the linear theory as a basis for estimating 
and f^ in which case they would not 'show a dependence on 6 . In 
such cases, the integrals in Eq. (33) become constants. When it can be 
further assumed that (t) ,oj is also independent of 6, the 

expression for takes the simple form 

Cjn-(6(t)) = A + BCn,^(»;6(t),o) (34) 

Despite differences in detail in the formation of the individual 
indicial responses, the representations of given in Eqs. (33) 

and (34) should find applications to other more general categories of 
wings and bodies. 

5.4 Experimental Evaluation of Stability Coefficients 

Experimental evaluation of requires a true simulation of the 

pitching motion illustrated in Fig. 6(b); attempts to simulate this 

motion have proven so difficult as to be impracticable. On the other 

hand, an experimental determination of the indicial pitching-moment 

response, and thus Cmg» from the motion illustrated in Fig. 6(a) 

requires instrumentation having such rapid response characteristics as 

to be impracticable. The interpretation of given previously might 

S 

provide a basis for an experimental determination of that coefficient 
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by means of appropriate energy measurements. However, difficulties 
remain in simulating the required motion, and in any case the interpre- 
tation does not suggest a basis for an analogous measurement of the 
more important coefficient . 

Rather than undertake these very difficult separate experiments, 
it can be argued that only one need be undertaken involving a combination 
of the terms and , since they usually appear together in the 

equations of motion. While this need not be true in general, it is 
certainly the case when the aircraft is assumed to follow an essentially 
rectilinear flight path, i.e., when q « a. This is easily seen from 
Eq. (26), which can be rewritten in the form 

+ [q(t) - a(t) ] -|- (t) ,oj (35) 

The last term vanishes identically when the flight path is precisely 
rectilinear, and so may be neglected when departures from a rectilinear 
flight path are small. Hence, only a single measurement is required of 
the term + YOrng » which is recognized as the damping- in-pitch 

coefficient. Clearly, a measurement of the coefficient that is in 
conformity with the way the term was derived requires a rectilinear 
flight path (q = a) and a motion involving small pitching oscillations 
in a about a fixed o equal to the instantaneous value of a under 
consideration. 
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6 . 


NONPLANAR MOTIONS 


The treatment of nonplanar motions raises the question of the role 
played by coupling in a nonlinear formulation. By coupling is meant the 
following: in the analysis of, for example, a combined pitching and yaw- 

ing motion, the linearity assumption (i.e., indicial responses independent 
of the past motion) allows the vector decomposition of the nonplanar motion 
into two orthogonal planar motions, the treatment of each planar motion as 
though the other were absent, and finally superposition of the separate 
results. In the nonlinear regime, where the indicial responses depend on 
all the variables of the past motion, the aerodynamic response to motion in 
one plane clearly will be influenced by the presence of the other motion. 
That influence is what is meant by coupling; the role assigned to it by the 
nonlinear formulation is studied in the sections to follow. 

For simplicity, flight-path properties, including flight velocity, 
are again assumed to remain constant throughout the motion, so that a 
general motion involves five independent variables. These are taken 

to be the components of the aircraft’s flight velocity and angular 

. 

velocity resolved either in the aerodynamic axis system (6, \p, X, q, r) 

/V 

or in the body axis system (a, 3, p„, q„, ^t,)* Having established the 
general integral form for the aerodynamic response to a motion involving 
two variables (Eq. 15), extending it to five variables is formal and 
requires no additional explanation. The general Integral forms are 
assumed to have been derived and simplified in accordance with the same 
simplification invoked in the analysis of the planar motions; namely, 
that the indicial responses within the integrals, initially functionals, 
have been replaced by functions dependent only on the magnitude of the 
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past motion Just prior to the origin of the steps. The analysis begins, 
then, with simplified integral forms analogous to the form applicable to 
planar motion, Eq. (18). 

6.1 Approximate Formulation in the Aerodynamic Axis System 

As noted above, the aerodynamic force and moment components 
resolved in the aerodynamic axis system are assumed to depend on the 
five variables 5, ij;, X, q, r. For example, the pitching-moment 
coefficient specified as a functional of the form 

C^(t) = G[6(5),ij^(5),A(5),q(?),r(5)] (36) 

The simplified integral form analogous to that of Eq. (18) thus contains 
five integrals, one for each variable. The indicial functions within 
the integrals depend on t - t and the magnitudes of the five variables 
just prior to the origin of the steps. Just as in the planar case, the 

integrals can be further simplified by reducing them to forms correct 

to the first order in frequency, assuming that 6 and if; may be large 
but that the angular rates X, q, r remain small. If the procedure 
illustrated for the planar case (Eqs. 19—27) is followed, an expansion 
of the integrals about X=0, q=0, r=0 yields, to first order in 

the rates, a sum of stability coefficients. The result is 

Cm(t) = C^(oo;6(t),iJ;(t)j + ^ C^* (oo;6 (t) 

. • 

+ T Cm6(«(t),Kt)) C„,^(6(t),i|,(t)) (37) 
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where, as before, the infinity symbol indicates steady flow. For brevity 
the zeros belonging to X, q, r have been omitted. Analogous expres- 
sions for and C and the axial, side and normal-force coefficients 

I n 

C„, are obtained by substituting these coefficients wherever C 

X Y Z ^ o jjj 

appears in Eq. (37). 

6.2 Simplification of the Formulation in the Aerodynamic Axis System 

Just as in the planar case (Eq. (35)), an additional simplification 
of Eq. (37) can be achieved by invoking the conditions of an almost 
rectilinear flight path. The conditions require, as before, q^^o, and, 
in addition, r w eA . Adding and subtracting the terms (c^2./V)Cjjj^ and 
e(XZ/V)Cui^ and rewriting Eq. (37) yields 

" C^(~;6(t),<|;(t)) +^C^^(6(t),ij;(t)) 

■•■■V jcmq(“;'S(t),<f'(t)j + YCm^(<S(t),i//(t))j 

+ (q - o) ^ ,i|j(t)j 

+ (r - eX) Y (t) ,ijj(t)j (38) 

The last two terms vanish identically when the flight path is precisely 
rectilinear, and so may be neglected when departures from a rectilinear 
flight path are small. The remaining terms are identified by comparing 
them with those obtained when the flight path is precisely rectilinear, 
where exactly, q = a, r = 6(j), X = y(j), and (J) is the oomng rate of the 
longitudinal axis around the flight velocity vector. The result is 
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^mq + YCmg 


( 39 ) 



+ ^C]j) 


Each of the terms in Eq. (38) is associated with a particular motion 
from which it may be evaluated. The term (t) ,i|;(t)j is the 

pitching-moment coefficient that would be measured in a steady planar 
motion with 6 and ij; at the fi:.ed inclinations 6 = const, ip = const. 
The term is a damping- in-roll coefficient that would be measured 

from small oscillations in ip about ip = const with 6 fixed at 
6 = const and <j) fixed at zero. The term is, as before, the 

damping- in-pitch coefficient for planar pitching oscillations about an 
axis normal to the plane of a, measured now, however, with both 6 
and ip at the fixed inclinations 6 = const, ip = const. The term C^* 

<t> 

• • 

is the rate of change with (pl/V, evaluated at of the pitching- 

moment coefficient that would be measured in a steady coning motion, 

6 = const, = const, cp = const. Thus, coning motion emerges as one 
of the characteristic motions required in constructing the response to 
an arbitrary nonplanar motion. 

• • 

In summary, with terms multiplied by (q - a) and (r - eX) neglected, 
the aerodynamic force and moment system takes the form 


C,^(t) - c^(. 



;6(t),t(t)) Cfc-(6(t),*(t)) 


= iX,Y,Z 
(l,m,n 


(40) 


Equation (40) suggests that, for nearly rectilinear flight paths, the 
forces and moments due to an arbitrary motion may be compounded of the 
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contributions from four simple motions: steady resultant angle of 

attack, oscillations in roll and pitch at constant resultant angle of 
attack, and coning at constant resultant angle of attack, all at a 
constant inclination of the body axes from the crossflow velocity 
vector. The motions are illustrated schematically in Fig. 10. 



Fig. 10. Characteristic motions in the aerodynamic axis system. 
Linear dependence on angular rates. 

6.3 Approximate Formulation in the Body Axis System 

The aerodynamic force and moment components resolved in the body 

axis system are assumed to depend on the five variables a, 3, q„, 

B B 

r_. For example, the pitching-moment coefficient C (t) is said to be 
n m 

a functional of the form 

C^(t) = H[a(5),3(5),Pg(5),q3(5),r^(?)] (41) 


The reductions of the integral form based on Eq, (41) parallel those 

of the preceding sections. Expanded about p„ = 0, q_ = 0, r_ = 0, the 

B B B 

/V 

integral form for yields, to first order in the rates 
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^ Cmp^(co;a(t) ,g(t)) 


q ^ r Z' 

+ ^ Cm (~;a(t),e(t)j +-|-Cm (~;S(t),6(t)) 


Cm^(»(t).^(t)] +-^ Cmg(a(t),3(t)) 


( 42 ) 


where the zeros belonging to p , q , r have been omitted. Analogous 

i3 is D 

expressions for and and the axial, side, and normal-force 

■^ ^ /\ 

coefficients C^, C^, are obtained by substituting these coefficients 
wherever appears in Eq. (42). 

That Eqs. (37) and (42) yield compatible forms may be verified by 

/s /\ 

transferring and to the aerodynamic axis system by the use of 

Eq. (9) and then replacing the variables a, B, p , q , r by variables 

D 

in the aerodynamic axis system through the use of Eqs. (4), (6), and (8). 

It will be found that each coefficient in C (t) and C (t) can be 

m n 

matched with a combination of coefficients in C (t) and C (t) having 

m n 

the same multiplying variable. The matches for the coefficients in 
C^(t) yield: 


Cjj^(“;6 ,i|^) = C^(«;a,3)cos ip - C^(~;a, 3)sin \p 


C = Cm («;a,3)cos ii» - C_ (~;a,3)sln ip 

m-i Pg 


pb 


Cma(“;'5,'J^) = Cm^_(“;oi,3)cos2 rp + Cn^_ (co;a,3)sin2 ijj 




B 


(43)' 


- [cnqB(“5“»^) + Cmj.B(“J“» &)] cos \p sin tp 


*Eq, (43) continued on following page 
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= Cm^g(»;S,3)cos2 ip - Cn^^(«>;aJ)sin^ ip 
+ jemq^(“;ot,§) - Cnrg(“;«.B)]cos ip sin ip 
Cia^(S,^) = ^^(a»3)cos2 4 , - Cng(a,e)sin2 ip 


+ [Cnig(ct,3) - Cn^(a,B)]cos \p sin ip 

Cm^(<S,'l^) = Cmp^(~;S,3)cos ^ - Cnp^(~;a,B)sin ^ 

+ 6Cm|(a,3)cos^ + 6Cn^(ct,3)sin2 \p 

- sj^Cm^Ca.e) + Cng(ot,3)]cos \p sin ij; 


( 43 ) 


The analogous matches for the coefficients in obtained 

from Eq. (43) by replacing with and with “fim^* 

matches for the damping-in-pitch coefficient and the side-moment 

coefficients Cn^ and are of particular interest. 

Equations (39) and (43) yield 


Cnip(*^>'J^) “ Qniq 


- + re.,i)cos2 * + - rC„jj8ln2 ♦ 


S, + (C«.3 - 'rC„g)J 


COS 4* aia ip (44) 
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Cnl(“;<S,4') = YCn^ + SCn^. 


= cos Tp + Cmp^ sin ip^ 


'1' + ^ 


+ 6 Cn CO 

+ (^nq^ + '*|] 


(45) 


Cn* - YCn,;, = 6{/C 


m, 


'iB 


YCni^)sin2 ip 


+ - YCng^cos2 ,p 


+ 


(e„,3 * vc„j) . - .e„:j 


COS \p sin ip] 


(46) 


Combining Eqs. (44) and (46) yields the following interesting equality: 

(47) 


Cn} - rc„j + 6C^ - + VC„:) + (c„^^ - yC^g) 


Cjn + YCrn^ I is recognized as being the planar damping-in-pitch 

qg “/ 

yN 

coefficient measured along y for small oscillations in a about 

j5 

a = const with 3 held fixed at 0 = const. Similarly, 
is the damping- in-yaw coefficient measured along z for small 
oscillations in 3 about 3 = const, with a held fixed at a = const. 
Thus, a measurement of would be equivalent to a measure of 

a combination of the three damping coefficients. The identity is shown 
schematically in Fig. 11. Equation (47) generalizes to the nonlinear 
case and to arbitrary bodies the relationship between - YCn^j 
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§ " 

(Cmq^+XCmA) + ( " ^'^<r 

Fig. 11. Schematic representation of the equality between 
- yCn^j/6 and the three damping coefficients. 

and the damping coefficients that was pointed out and verified for 
bodies of revolution in the linear case in Ref. 5. It is noted that 
(Cn^ - equals 6 ^Cmq^ + a = 0 and equals 

6 - yCng) when 3 = 0. Under conditions where a linear formulation 

of the moment system can be assumed to hold (e.g,, when a 0, 0-^0) 
it is consistent to assume that the couplings between motions in a 
and 3 will be negligibly small. Under these conditions, the measure- 
ment of at a = Q and again at 3=0 is all that is 

required to yield measures of the damping coefficients characteristic 
of the two uncoupled modes. 

6.4 Simplification of the Formulation in the Body Axis System 

When the assumption of a nearly rectilinear flight path is justi- 
fied, Eq. (42) can be simplified to yield a form analogous to that of 
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Eq. (38) in the aerodynamic axis system. The approximate expressions 
for qg and r^, consistent with q a, r «s eX, are 


a . 8 

q., « — + P., — 
Y Y 






Substituting In Eq. (42) gives 


in 


C (~;a,B) 

in 


1 Vr 

Y V L' 


YCmn (”;a.e) + 


(“’;«, 8) 

‘iB 


+ YCm^(a,6)j -YCmg(a,8)j 


(49) 


Again, each of the terms in Eq. (49) is associated with a particular 
motion from which it can be evaluated. The first term is the pitching- 
moment coefficient along y_ that would be measured in steady planar 

a 

^ ^ ^ /V 

motion with a and 3 at the fixed inclinations a = const, 3 = const. 
The combination of terms multiplied by p Z/Vy can be shown (from 

u 

Eq. (45)) to be the rate of change with ^l/Vy evaluated at (|) = 0, of 
the pitching-moment coefficient along y_ that would be measured in a 

a 

^ /s . 

steady coning motion a =* const, 3 - const, <j) = const. The third term 
has already been defined as the damping- in-pitch coefficient measured 
along y for small oscillations in a about a = const with 3 held 

D 

fixed at 3 = const. The term ^ cross-coupling term 

resulting from the oscillatory motion in the ^ plane required for the 
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evaluation of the damping- in-yaw coefficient. This term and the 
analogous term in /^nq ■*" ^^n^j nonlinear interaction 

terms that are normally excluded in the classical treatment and are 
missed by attempts to generalize from linear formulations based on the 
principle of superposition. 

In summary, for nearly rectilinear flight paths, the aerodynamic 
force and moment system in body axes take the form 


1 


7 ^ - YCki(5.6)j ; k - 




m ^ /s. 


(50) 


In the body axis system, the four characteristic motions are steady 
angle of attack and sideslip, coning at constant angle of attack and 
sideslip, and the oscillations in pitch and in yaw at constant angles 
of attack and sideslip. The oscillations-in-roll motion that was 
required in the aerodynamic axis system is, in effect, incorporated in 
the oscillations in pitch and yaw motions in the body axis system. An 

^ A*. 

oscillation in a with 3 held fixed, for example, will be seen to 
involve an oscillation in The four motions are illustrated 

schematically in Fig. 12. 





STEADY ANGLE OF ATTACK YAW OSCILLATIONS PITCH OSCILLATIONS 
AND SIDESLIP 

Fig. 12. Characteristic motions in the body axis system. 
Linear dependence on angular rates. 
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6. 5 Body of Revolution 

Symmetry considerations in the body axis system can be used to 
further specify the form of the aerodynamic coefficients. The arguments 
will be carried out for the body of revolution, since they can be carried 
farthest in this case and will indicate the form of the interaction 
terms . 

Consider an arbitrary motion in a and g with p_ = 0. The most 

B 

general form that can be written for the pitching moment that is of 
third order in the angles and first order in the rates is 

^m^^^ ^0 ^ a2a^'+ a3U^ + a4§ + asg^ + + ayag + a0a^g + agag^ 

V (ho + bxa + b20t^ + b3g + b 40 ^ + b5ag) 

+ -^-^(cQ + cia + 02^^ + 036 + + csa 3 ) ( 51 ) 

The requirement that C be an odd function of a and an even function 

m 

of g reduces Eq. (51) to 

• • 

Cuj(t) = a(ai + a 3«2 + ag3^) (bg + b 2 a 2 + b4§^) +^-^(c5ag) 

(52) 

The requirement that C be equal to -C when a and g are inter- 

n m 

changed gives for C (t) 

n 

• « 

-C^(t) = 3(ai + a3B^ + aga^) + “ (bg + b20^ + b4a^) (cgag) 

(53) 
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Finally, the requirement that the total moment in the crossflow plane 
be normal to the direction of 6 when the motion is planar (ij) = 0 ) 
gives 


as - ag 
b4 = b2 “ C 5 



so that 

• • 

C^(t) = a(ax + 336 ^) + Y (bo + i>2^^ ~ 056 ^) + — ^ (cgag) 

• • 

-C^(t) = g(ai + a362) + i ^ (bg + ba^^ - cgS^) +^f- (cgSe) (55) 


Matching coefficients in Eq. (50) yields 


"m 

(a>;a,8) = 

a(ai + asd^) 



+ = 

(ho + 


058^) 

Cm^ 

- = 

yv As 

-Cgae 



c 

n 

(~;a,6) = 

-6(ai 

+ as62) 


A • 

- YCng = 

(bo + 

- 

^2\ 

csa^) 

y\ 

Sb 

^ ♦ 

+ YCn- = 

/\ 

-CgaS 




(56a) 


\ 


(56b) 


For the body of revolution, the interaction terms are equal and their 
form to second order is a product of a and Equation (56) can be 
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used to find the analogous terms in the aerodynamic axis system. 
Equations (43) — (47) yield 


Cjjj(“;<S,'/') = 6(ai + 

Cn^(6,'j') = (bo + b25^) 

^m^ “ = 0 

Cj^(“;«S,'l') = 0 

Cn^ («.’*') = 0 

Cnj, ~ = 6(bo + - 056 ^) 


(57a) 


(57b) 


As required, the coefficients in the aerodynamic axis system are 
functions only of 6. Note that C„* equals /Cm + yC™M only 

“XT \ “^06/ 

when § = 0 and equals only when a = 0. That this must 

be true may be verified by inspection of Fig. 11. The relationship 
between and - YCn^^/<5 is of particular interest. Equa- 

tion (57) indicates that to first order in 6 (i.e., small 6 where a 
linearized theory can be expected to hold), the two terms must be equal; 
but this requirement does not hold for larger values of 6 where terms 
of 0(6^) must be retained. The consequences of the breakdown of the 
equality have been discussed several times (cf«, for example, Refs. 6 
and 14). Briefly, breakdown of the equality is proof that the linear 
stability derivative formulation cannot be extrapolated into the 
nonlinear regime simply by allowing the coefficients that appear as 
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I 


identical constants in the linear formulation to become identical 
functions of 6. A fomulation incorporating this erroneous assumption 
can be shown to yield misleading results when used in a program for 
extracting nonlinear aerodynamic coefficients from free-f light data [6]. 

The significant aspects of the relationship between 3.nd 

(Cn^ - have been confirmed by the results of two recent 

experiments. The first involved spinning and coning motions of a slender 
cone at a supersonic speed [ 5 ], The results, obtained with an apparatus 
that allows investigation of separate or combined spinning and coning 
motions, are important in two respects: first, experiments with a 

purely spinning cone at constant 6 (i.e., the classical Magnus experi- 
ment 6 = const, ^ = const, = 0, which, for bodies of revolution, can 
be used to replace the oscillations-in-roll experiment required to 
measure for nonaxisymmetric bodies) failed to reveal a measurable 

Magnus moment coefficient and second, the results confirmed that, 

with negligibly small, the equality that should exist for small 

6 is an equality between the side-moment coefficient due to coning 
Cj^^/6 and , The main results are shown in Fig. 13 . Measured 

values of Cn^ are seen to be in excellent agreement at small 6 with 
the straight-line variation obtained from a linear potential 

theory [ 15 ]. The breakdown of the equality at larger 6 has been 
confirmed by the results of a second experiment carried out by 
Iyengar [ 16 ]. This was a careful oscillations-in-pitch experiment for 

over a large range of 6 with both model and wind-tunnel conditions 
the same as those used in the coning experiment. The results, also 
shown in Fig. 13 , clearly confirm that the equality does not hold at 
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s 


Fig. 13. Variation of with 6 ; 

Mach number - 2.0 y I /t = 0.61. 

eg 


large 6 . Equation (57) shows that the source of the inequality between 

and - yCn^j /6 lies in the term C 5 , the coefficient of the 

cross-coupling terms ~ ^^mg j ( CnqB ’Y^n^) “ The experimental 

results indicate that |c 5 | > |b 2 |> where b 2 is the coefficient of the 

nonlinear contribution to . That is, the cross-coupling terms are 

at least of the same order of magnitude as the nonlinear part of the 

principal damping coefficient, and so must be retained for consistency 

whenever nonlinearities are admitted into the principal terms of the 

aerodynamic formulation. Why the cross-coupling terms must exist when 

a and 3 become sufficiently large is indicated below. 

Consider a body of revolution in uniform sinking and sideslipping 

motion a = const, 3 = const, with p„, q^,, r = 0. The motion is 

B 

illustrated in Fig. 14. The resultant normal-force coefficient 

^ /V 

Co(“;ot, 3 ) initially is alined with the resultant angle-of-attack vector. 
K 

^ A. 

At time zero, the body is given a step change in a with 3 held fixed. 
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TIME 0 


TIME t 


TIME — CO 



Cr (coja, / 3 ) 


Fig. 14. Maneuver showing origin of cross-coupling term 

* n 


This is the experiment that would yield the indicial normal-force 
response C 2 ^(t;a, 3 ) to which the term Cz^ is related by an integral. 
As shown in Fig. 14, subsequent to the change in a the resultant 
normal-force coefficient Cj^(t;a + Aa,3), in addition to growing, must 
rotate in order that it be alined with the resultant angle of attack 
again after a steady state has been established. Clearly, a step change 
in a will induce a time-dependent side force in the interval during 
which the resultant normal force rotates to its new position. The 
induced time-dependent side force is related by an integral to the 
cross-coupling term Cy* that appears in the expression for C^Ct). 

/V 

Such a term must exist whenever a maneuver in a occurs in the presence 
of a nonzero 0 . 
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6.6 Aircraft Spin Motions 

The emergence of coning motion as a characteristic motion in both 
the aerodynamic and body axis systems emphasizes its importance in the 
nonlinear formulation. With p_ = 0, coning motion is equivalently the 

D 

simultaneous periodic variation of two orthogonal planar motions; in 
thus enabling two planar motions to interact, coning motion can be said 
to characterize the coupling problem. 

Further evidence of its important is the obvious similarity between 
coning motion and the steady spin of an aircraft, suggesting that a 
moment formulation based on either Eq. (40) or Eq. (50) will properly 
describe the aerodynamic responses of spinning aircraft. This view is 
encouraged by the success achieved in two attempts to reproduce aircraft 
or model spin motions by calculations based on aerodynamic formulations 
bearing a similarity to those proposed here. In the first [17], the 
actual spins of an F-lOO aircraft were reproduced by calculations based 
on an aerodynamic formulation that called principally for wind-tunnel 
measurements of the conventional static forces and moments. In the 
second [18], the spins of a delta-wing model in a spin tunnel were 
reproduced by calculations based on an aerodynamic formulation that 
called principally for wind-tunnel measurements of the forces and moments 
on a model in coning motion. The formulations based on Eqs. (40) and 
(50) in effect include these motions within the four characteristic 
motions whose force and moment contributions are required to build up 
the response to an arbitrary motion. It is known, however, that in the 
establishment of a spin the large asymmetric regions of separated flow 
on the wings of the aircraft may cause the aerodynamic responses to be 
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nonlinear functions of the spin rate, even at low spin rates. This 
contradicts the assumption underlying the development of Eqs. (40) and 
(50) that the aerodynamic responses be linear functions of the rates. 

The authors have shown recently [9] how the formulations can be general- 
ized to allow a nonlinear dependence on the coning rate, which should 
make them more fully applicable to the analysis of spin motions. The 
restriction imposed previously on the degree of cognizance of the past 
motion allowed the indicial response remains in force, however, so that 
aerodynamic hysteresis effects still cannot be acknowledged within the 
scope of the otherwise extended formulations. The main results are 
indicated below. 


6.6.1 Aerodynamio Axis System: The generalized formulation paral- 
leling that of Eq. (40) (again assuming a nearly rectilinear flight path) 
takes the form 






_ X,Y,Z 




(58) 


That Eq. (58) reverts to Eq. (40) when the linearity condition on coning 
rate is reimposed can be seen by expanding the terms in Eq. (58) to first 
order about X/y = 0. The first term in Eq. (58) is the coefficient that 

would be measured in a steady coning motion 6 = const, ip = const, 

• • 

(p = X/Y = const. As before, the term is the damping- in-pitch 

coefficient that would be measured from small oscillations in a about 
a = const with fixed at \jj - const but now, in addition, in the 
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presence of a steady coning motion I = A/y = const. Similarly, is 

the damping- in— roll coefficient that would be measured from small oscil- 
lations in about ij; = const with 6 fixed at 6 = const and in the 
presence of a steady coning motion ^ = A/y = const. The indicated 
functional dependence on 6, A/y must be interpreted as follows: 
for flight with given values of 6, i|;. A, q, r at a particular instant, 
the aerodynamic coefficients that are to be associated with that instant 
are those evaluated around a coning motion having constant values of 6 
and ijj equal to the instantaneous flight values and a constant value 
of coning rate equal to the instantaneous value of A/y. 

Thus, the four contributions required in Eq. (40) to build up the 
response to an arbitrary motion reduce to three when a nonlinear 
dependence on coning rate is admitted. This is because the first term 
in Eq. (58) C^(<»; 6 ,ij;, A/y) is the general term which, it now appears, 

replaces two terms in Eq. (40) representing the expansion of 
Ck(°°; 5 A/y) around A/y = 0 to first order in A/y. The more important 
change, however, at least from the experimental standpoint, is that 
retaining a nonlinear dependence on coning rate requires for consistency 
that the oscillatory experiments be carried out in the presence of coning 
motion. The three motions are illustrated schematically in Fig. 15. 



AND CONING AND CONING 


Fig. 15. Characteristic motions in the aerodynamic axis system. 
Nonlinear dependence on coning rate. 
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6.6.2 Body Axis System: The generalized formulation paralleling 


that of Eq. (50) takes the form 


Ck(t) = Cj^(”;a,6,Pg/Y) + 


lk\i 

Y V [■ 


Ckn (“;a»3,P„/Y) + 

qs ^ 


YCk^(ct,e.Pg/Y)j 


« 

- - YCkg(S,a.Pg/Y)J ; k = (59) 


Equation (59) reverts to Eq. (50) upon expanding the terms in Eq. (59) 
to first order about p„/y = 0. The first term in Eq. (59) is the 

D 

coefficient that would be measured in a steady coning motion 
^ ~ Ptj/y ” const with a and S at the fixed inclinations a = const, 

5 = const. The second term is the damping-in--pitch coefficient that 
would be measured from small oscillations in a about a = const with 

6 fixed at 0 = const and in the presence of a steady coning motion 

• 

^ = Pt,/y ” const. The third term results from small oscillations in 

0 about B = const with a fixed at a = const and in the presence 

of a steady coning motion 0 = P-d/y* It should be noted that the coning 

rate j) = p^/y on which the terms depend in Eq. (59) is not equal in 

magnitude to the coning rate ^ = A/y in Eq. (58). The rates differ 

by ^/y. Thus, in the body axis system, the appropriate constant value 

of the coning rate for the aerodynamic coefficients that are to be 

associated with an instantaneous flight condition is that formed from 

the instantaneous value of p„/y. The three motions required in the 

B 

body axis system are illustrated schematically in Fig. 16. 
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Fig. 16. Characteristic motions in the body axis system. 

Nonlinear dependence on coning rate. 


6.6.3 Spin Radius: Equations (58) and (59) were derived on the 

assumption of a nearly rectilinear flight path, which would apparently 
restrict their application to spin motions having essentially zero spin 
radius. It can be shown, however, that the results will apply as well 
to spin motions having constant spin radius. Motions having a constant 
spin radius can be characterized by the existence of a point other than 
the mass-center about which the body rotates. This point, which lies on 
the body x^ axis, is itself in essentially rectilinear motion. The 
existence of such a point usually will guarantee fulfillment of the 
conditions under which Eqs. (58) and (59) apply, namely, q - a 0, 
r - eX « 0. The principal restriction is that | qx/V | , | rx/V | << 1, 
where x is the distance along x between the mass-center and the 
point in nearly rectilinear motion. Variations in spin radius x sin a 
also can be tolerated under the additional restriction |x/v| << 1. 
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6. 7 Requirements for Experiments 

Within the assumption of a linear dependence of the moment on coning 
rate, the analysis suggests that the moment contributions resulting from 
four characteristic motions are required to completely specify the 
nonlinear moment system for arbitrary motions about nearly rectilinear 
flight paths (Figs. 10 and 12). For wind-tunnel tests in the aerodynamic 
axis system, two kinds of apparatus would appear to be necessary: 

(1) A coning and spinning apparatus similar to the one described 
in Ref. 5. (For nonaxisymmetric bodies, the spin motor that 
reproduced the constant spin rate ^ would have to be replaced 
by a device reproducing small oscillations in about 

= const.) Such an apparatus should be capable of measuring 
the moment contributions due to steady resultant angle of 
attack, coning at constant resultant angle of attack, and 
oscillations in roll at constant resultant angle of attack. 

(2) An oscillations-in-pitch apparatus for measuring the moments 
due to small oscillations in a about a fixed a with the 
axis of rotation oriented normal to the a plane. It is 
emphasized that for nonaxisymmetric bodies, the oscillations- 
in-pitch apparatus must be capable of measuring not only the 
pitching moment but also the induced side and rolling moments. 

Wind-tunnel tests in the body axis system require the same coning 
experiment as in the aerodynamic axis system and, in addition, separate 
oscillations-in-pitch and oscillations-in-yaw experiments. The 
oscillations-in-pitch device must be capable of simultaneously measuring 

/V 

all three moment components due to small oscillations in a about a 
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/s. 

fixed a with 3 held fixed. The same device can be used in the 
oscillations-in-yaw experiment, in which the roles of a and 3 are 
reversed. Hence, experiments carried out in the body axis system, 
requiring only a single oscillatory device, would appear to have an 
advantage over those in the aerodynamic axis system from the standpoint 
of economics. It is essential, however, that the oscillatory experiments 
in the body axis system enable the measurement of the nonlinear inter- 
action terms even for bodies of revotution. Recently, the successful 
development of a device capable of measuring these interaction terms has 
been reported [19]. 

In programs designed to extract the nonlinear aerodynamic 
coefficients from free-flight data, of course it is recommended that a 
form for the aerodynamic moment system based on Eq. (40) or Eq. (50) be 
incorporated. It has already been noted [6] that procedures not allowing 
for the presence of interaction terms in the representation of the moment 
system can assign erroneous weights to the remaining terms. 

Eliminating the assumption of a linear dependence of the moment on 
coning rate reduces the number of characteristic motions required from 
four to three: in either axis system, a coning motion and two 

oscillatory motions in the presence of coning (Figs. 15 and 16). Experi- 
ments designed to reproduce the motions in the wind tunnel again require 
the coning apparatus and the types of oscillatory devices just described. 
The significant additional requirement that each of the oscillatory 
experiments be carried out in the presence of coning means, of course, 
that now the oscillatory devices must be incorporated in the coning 
apparatus. These obviously difficult experiments, involving oscillatory 
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and coning motions in combination, are required only where the moment 
contribution due to steady coning shows a significant nonlinear depend- 
ence on coning rate. Otherwise, the experiments may be conducted 
separately as described above. Useful surveys of the needs and 
capabilities for carrying out experiments in the wind tunnel involving 
oscillatory, rotary, and combined motions at high angles of attack are 
available in Refs. 20 and 21. 

7. CONCLUDING REMARKS 

A study aimed at extending the notions of aerodynamic indicial 
functions and superposition integrals into the nonlinear regime has shown 
that replacing the indicial functions within the integrals by functionals, 
themselves dependent on the past motion, achieves the desired objective. 

A simple specification of the degree to which the reformulated indicial 
responses are cognizant of the past motion led to practicable nonlinear 
generalizations of the linear superposition and stability derivative 
formulations of the aerodynamic response to arbitrary motions. Applied 
to arbitrary nonplanar motions, the generalization yielded a form for 
the aerodynamic response built up of the contributions from a limited 
number of well-defined characteristic motions, in principle reproducible 
in the wind tunnel. Further generalizations, enabling the acknowledg- 
ment of more general categories of aerodynamic phenomena within the 
scope of the formulation, are possible. The characteristic feature of 
hysteresis, multivalued aerodynamic responses, could be acknowledged, 
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for example, by extending the degree to which the indicial responses are 
cognizant of the past motion one step beyond that specified in this 
study. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., 94035, May 1975 
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